Twin stimulated amplification of light and matter waves in an atom— photon pair laser 
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We consider ultracold atoms in a far detuned optical lattice orientated across a high-Q optical res- 
onator. Applying an external driving laser to the atoms, which is red detuned from the cavity mode 
by one vibrational quantum, induces cavity-enhanced sideband cooling. For a dense and cold enough 
atomic ensemble we predict an oscillation threshold for optical Raman sideband lasing concurrent 
with coherent matter-wave amplification. Above this oscillation threshold photons and atoms in 
the lowest band are dominantly created pairwise via stimulated emission with a strong suppression 
of competing spontaneous processes. In close analogy to a nondegenerate parametric oscillator we 
find sub-Poissonian photon statistics and almost perfect nonclassical atom-photon number correla- 
tions. Injecting atoms in higher vibrational bands via tunneling or incoherent scattering then leads 
to continuous, simultaneous generation of a coherent atom beam and laser light with nonclassical 
atom-field correlations. 



I. INTRODUCTION 

Already a century ago the striking analogy between 
atomic Schrodinger waves and electromagnetic waves has 
started a new field in physics studying wave properties 
of massive particles [T|. Nevertheless, only the significant 
advances in the cooling and trapping techniques of neu- 
tral atoms in the last decades 0] made extensive and well 
controlled experimental investigations of macroscopic co- 
herent matter waves possible 

A central and ongoing quest in this field is the realiza- 
tion of the atomic counterpart of an optical laser generat- 
ing a coherent intense matter-wave beam. Although first 
proposals to implement such an atom laser suggested the 
use of refined optical cooling techniques within optical 
traps 0, S , it was the creation of Bose-Einstein con- 
densates via evaporative cooling which paved the way to 
coherent atomic beam sources by output coupling from 
such a condensate 0, i, i, 13 • Plenty of follow -up exper- 
iments based on this principle have been performed since, 
including the study of the spatial I 111 as well as tempo- 
ral [l3| coherence, the divergence and the counting 
statistics [l^ of the produced atomic beams. 

Despite some improvements a quasi continuous-wave 
(cw) implementation of an atomic beam out-coupled 
from a pre-existing Bose-Einstein condensate (BEC) is 
still hampered by the finite number of atoms initially 
present in the condensate. This limits the pulse duration 
and the phase and intensity stability of the matter-wave 
beam. The central challenge to overcome here is con- 
tinuous reloading of atoms. Naturally a pump source is 
needed for a cw operation of a laser. As atoms cannot 
be created out of nothing, one cannot just inject energy 
itself but the only way to do this for an atom laser is to 
feed new atoms into the BEC mode [HI, [1^1 ■ 

In a straightforward but technologically challenging 
approach, one has tried to convert the time sequence of 
BEC generation into a linear array of spatially separated 
cooling stages to operate a BEC generator in the contin- 
uous regime [l^l- These efforts already produced a very 



cold and partially coherent atomic beam, but to achieve 
full atom lasing, the setup still needs further improve- 
ments. 

Theoretically, several alternative proposals for a con- 
stantly operating atom laser by continuous atomic out- 
coupling from a BEC were put forward including continu- 
ous injection of cold atoms and ongoing evaporative cool- 
ing [l3i0,ll3|- The re-filling of condensate atoms is then 
achieved by a coUisional transfer of atoms into the BEC 
mode out of a large reservoir. While the corresponding 
microscopic implementation of the injection mechanism 
may be manifold, the temperature of the newly injected 
atoms typically has to be close to the critical value for 
condensation. When using near-resonant lasers for the 
required optical pre-cooling during the evaporative cool- 
ing process, one again recovers the problem of reabsorp- 
tion heating [2l| as only high phase-space densities can 
assure the required atomic collision rates for fast enough 
replenishing. These problems are then closely related to 
the difficulties appearing in the early proposals for atom 
lasing via all optical cooling d, 0, [5^ . 

In this paper we study a different scenario for contin- 
uous optical replenishing the BEC mode at much lower 
density by help of an optical resonator to tailor the op- 
tical atomic cooling process Such a cavity helps to 
suppress spontaneous emission and thus indirectly pre- 
vents reabsorption. Naturally cavity cooling in its stan- 
dard form was already considered as a tool to achieve a 
BEC by all optical means [23|. In general the final tem- 
perature is only limited by the cavity linewidth, and pho- 
tons emitted through the mirrors cannot be reabsorbed. 
Hence, this technology should be in principle sufficient to 
achieve degeneracy, if a good enough cavity can be build 
around a dense atomic cloud. Such a very narrow cavity 
line requires however a challenging technological setup, 
and the timescales involved in the dynamics get rather 
slow. 

Here we study an improved and conceptually simpler 
variant of a cavity-based setup for an atom laser, which 
is designed to work under different and technologically 



2 



less stringent conditions. The basic idea is to implement 
a combination of an atom laser and a photon laser op- 
erating simultaneously on a vibrational anti-Stokes Ra- 
man transition. Here the central physical mechanism to 
suppress the unwanted transitions and photon emissions 
brought about by the cavity is not strong atom-field cou- 
pling but stimulated enhancement of the desired transi- 
tion. This allows to suppress reabsorption as well, but 
still implies fast transition times. Note that a related 
setup was suggested to enhance sideband cooling in ion 
traps [2^, [2^ . Instead of a single trapped ion, here we 
consider a large ensemble of ultracold atoms in a far de- 
tuned optical lattice and the case of a good cavity. 

The occupation of the lowest lattice states is then in- 
creased by a stimulated transfer from a higher vibrational 
excitation (or a band in a lattice). Actually, for a single 
trapped atom, coherent coupling of atomic quantum mo- 
tion and a field mode in cavities via stimulated Raman 
transitions has been suggested and investigated already 
in earlier work [23] ■ In an anti-Stokes transition, even 
the vacuum field can stimulate the transfer of an atom 
into a vibrational lower state, simultaneous with induced 
photon emission into a cavity mode. When the photon 
leaves the resonator no backward process is possible. 

For a larger ensemble of N atoms this single photon 
emission process is collectively enhanced by a factor N 
with respect to a spontaneous transition [28[. Entering 
the regime where the Raman gain on the sideband transi- 
tion exceeds the cavity losses, more than a single photon 
is present and lasing via stimulated emission will start. In 
this case one can make use of double Rose enhancement 
through macroscopic occupation of the atomic mode and 
the cavity mode. Obviously the total transition rate is 
then proportional to the total photon emission through 
the cavity. 

Including replenishing of atoms to higher vibrational 
lattice levels and out-coupling of atoms from the low- 
est band should then lead to continuous lasing and atom 
beam emission. The goal of the present paper is to de- 
velop a consistent theoretical description of the essential 
physics of this pairwise stimulated enhancement and to 
study the nature of the generated output beams of light 
and atoms in more detail as already presented in our 
short letter Due to the complexity of a practical 

implementation we have to refrain to a series of simpli- 
fying assumptions to derive a solvable model as outlined 
below. 

From a different perspective the present paper can also 
be viewed as a continuation of the work on cavity-cooling 
of atoms induced by stimulated emission [2^, [s^, [3l| . 
There we investigated the light forces in a laser system 
composed of a single inverted atom [2^, [s^] or even an 
inverted atomic ensemble ^31j in a high-Q cavity. We 
found that the optical gain through stimulated emission 
in the system can be tied to motional cooling and atomic 
self-trapping. In this system several effects including las- 
ing, trapping, and cooling work together in an appropri- 
ate way such that atomic temperatures well below the 




FIG. 1: (a) Illustration of the cavity-lattice configuration 
forming the atom-photon pair laser. Both light and mat- 
ter fields are stimulatingly amplified by Raman-gain involving 
blue motional sidebands. 

(b) Illustration of the atomic motional dynamics including 
injection {^), cooling (7;), and out-coupling as well as 

losses due to heating {ntc, Hcb, i^c)- 

limit of passive cavity-cooling are possible (sol . [sH . Nev- 
ertheless the system stays far away from degeneracy as 
the atomic inversion is accompanied by noise and spon- 
taneous emission. Here by help of an external potential 
and employing a vibrational Raman transition, we can 
relax the condition of internally inverted atoms. In ad- 
dition, the external lattice provides for trapping so that 
the light force induced by the laser light are not needed 
for trapping, and thus significantly less dipole heating is 
present. 

It is well established theoretically and practically that 
sideband cooling of a single atom (ion) efficiently reaches 
the ground state. Theoretical predictions suggest that 
this behavior can be further improved utilizing a cavity 
[25l . [26l . [3^ . While in the free-space case generalization 
to many particles is a hard task, collective enhancement 
should even be helpful in the cavity-enhanced case [28j . 

The central remaining problem still is reloading of new 
atoms without disturbing the Raman gain process. There 
seem to be several solutions to this process. As a first 
suggestion the lasers applied for creating a lattice or a 
trap within the cavity can also be tailored to guide ex- 
tra atoms from a pre-cooled reservoir into the active las- 
ing zone as experimentally demonstrated [H, [H, [H, . 
This can happen in quite different forms. When lasers are 
applied to form, e.g., a quasi ID optical lattice potential 
along a standing wave laser beam, the guided atoms can 
simply enter and leave the cooling region on the basis 
of tunneling In another possibility the laser directly 
involved in the final cooling stage to reach degeneracy 
by pumping atoms from a cold reservoir into the lasing 
modes could be used [s^l . This even avoids the need 
of spatial separation of the pre-cooling and the lasing 
stages. 

Here, as sketched in Fig. [Ija), we consider a quite 
generic scenario where several of the above possibilities 
are included. Atom loading from a reservoir into the trap 
potential is simply modeled by effective rates, which can 
be coherent as for tunneling or incoherent via scattering. 
Our work is organized as follows: In Sec. [Til we give an 
overview of our atom-photon laser scheme. In Sec. IIIIl 
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we model it by an incoherent transfer of atoms into the 
trap as it could be realized, e.g., from a thermal source 
[l8l | or by an incoming cold atomic beam On the 

other hand, we use in Sec. IIVI a model based on atomic 
tunneling through the lattice which we already used in 
[2^. In Sec. mil we use a rate-equation approach for the 
excitation probabilities of the pair laser and study mean 
number and number statistics in the steady state. An ex- 
tended model description based on stochastic differential 
equations is derived in Sec. IIVI which allows us to inves- 
tigate the noise properties and the spectra of the output 
beams in Sec. |Vl 



II. MODEL 

Let us now present a mathematical model for a system 
as shown in Fig. [l] which exhibits the mechanisms un- 
derlying an atom-photon twin laser in a rather generic 
form, omitting as much specific experimental details as 
possible. The basic components are a large ensemble of 
cold two-level atoms placed in an external optical trap 
which is intersected by the field of an optical resonator 
with high quality. We consider the simple situation where 
the atomic trap consists of an optical lattice across the 
resonator. At each lattice site we assume at least two 
well separated vibrational energy levels. The manifolds 
of eigenstates of the lattice sites within the cavity are 
identified with the atomic source and lasing modes. In 
the simplest case of a long- wavelength lattice, only a sin- 
gle lattice site might be located within the resonator 
mode volume. Analogous to cavity enhanced sideband 
cooling [25| . we then tune an external pump laser and 
the cavity mode frequency to Raman resonance between 
these two trap states, which resonantly couples the cav- 
ity photons and the atomic (quantized) motional states 
|27| . We choose the frequencies in a way that an atom is 
transferred from the upper (source) to the lower (ground) 
vibrational level whenever a photon blue detuned to the 
pump is generated in the cavity mode. 

For large enough detuning from the upper Raman level, 
the only significant dissipation channel in this system is 
cavity decay. In this limit each photon emitted through 
the cavity mirrors increases the occupation number in the 
lower atomic state by one. The process in its ideal form 
only stops when no more atoms are present in the upper 
level |28(|. The central important point to note here is 
that the Raman coupling and thus the rate for the pair 
creation of photons and ground-state atoms is collectively 
enhanced with growing cavity-photon number and num- 
ber of atoms already present in the ground level. Hence 
even if it is weak for a single atom, for increasing num- 
bers of atoms and photons this desired process will more 
and more dominate other, unwanted transitions in the 
system. 

Note that in a realistic three dimensional configuration 
one probably will have several lattice wells within the 
cavity mode. Nevertheless, as long as collisions or atomic 



on-site interactions are weak, we can formally combine 
them into a single effectively trapping state or mini band. 
Hence we effectively replace the lattice bands simply by 
two effective levels (mini bands) describing the atomic 
source and laser modes 23]. 

The free Hamiltonian is just represented by three os- 
cillators in the form 

+ Eob^b + Eic^c. (1) 

The operators a, 6, and c represent the cavity photons 
and the two atomic modes describing the atom laser and 
atom source states, respectively; Wa denotes the bare cav- 
ity resonance frequency; Eq and Ei are the energies of 
the atomic bands. In a first approximation the Raman 
transition can then be described by an effective nonlinear 
three-mode coupling Hamiltonian of the form [ssj 

= ih-q {a)h^c- ahc^) (2) 

where 77 is the interaction strength. Note that with re- 
spect to our later goal of atom lasing, we have already 
dropped the non-resonant heating terms. For a single 
particle this is only valid for sufficient trap depth [2^ . A 
more realistic approach has to contain off-resonant tran- 
sitions as, e.g., discussed in Ref. [28] . 

In order to allow for steady-state operation, we have to 
include losses from the laser mode. Possible such coher- 
ent output-coupling processes include, e.g., atomic tun- 
neling [J] or the application of a Raman transition for 
an atomic state change into an untrapped state 113] . 
Mathematically, independent of the specific implementa- 
tion, the out-coupling mechanism is treated via a linear 
loss term in the master equation in Lindblad form anal- 
ogous to photon loss. The master equation including 
cavity decay reads thus 

p = -i/h [Hint,p] + Cp, (3) 
Cp = Ka{ [a, P + [ap, a'^] ) 

+ Kb{[b,pb^ + [bp,b^), (4) 

where Kq and Kb denote the cavity (amplitude) decay and 
atomic out-coupling strength, respectively. 

In addition to these terms we will also include spon- 
taneous Raman transitions from the upper to the lower 
atomic states without generating a cavity photon occur- 
ring at rate Kcbi which will be important only below 
threshold but are included for completeness here. The 
corresponding term Cct in the master equation reads 

Ceb^ n,b{[b^c,pch] + [bhp,c%]). (5) 

Note that losses induced by the Raman pump laser on the 
atomic ground and upper states can be simply included 
into effective corresponding loss rates («{,, Kc) and will be 
small for a far detuned Raman transition. 

A central and critical point in the model is of course 
the refilling of the atoms in the source mode, which is 
required for a steady laser operation. As the size of the 
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full lattice is typically orders of magnitudes larger than 
the cavity volume, the trapped atoms outside the mode 
already form a useful reservoir for this purpose. In princi- 
ple this reservoir might be even periodically refilled from 
a nearby molasses operated far detuned from the Raman 
pump laser. A more decisive challenge is to design a 
mechanism how these reservoir atoms then can be con- 
tinuously transferred into the interaction region within 
the cavity. While various methods for this are conceiv- 
able in practice, at the end they can be modeled by only 
two different generic types of pumping, namely via coher- 
ent or incoherent transfer of atoms to the upper source 
mode. While the first corresponds to a tunneling type 
injection from filled sites outside the resonator or a de- 
terministic injection by actively moving the filled lattice 
through the resonator [34] , the second describes refilling 
by scattering between atoms or optical pumping from 
higher lying states. 

Before discussing these two possibilities and their con- 
sequences in sections IIIII and IIVI in detail, we will first 
give some qualitative insight into the system via a di- 
rect numerical solution of our model, i.e., Eq. ([3]). This 
is only possible for very low atom and photon numbers, 
but nevertheless some general features of the dynamics 
can be already nicely seen here. In Fig. [D^a) we plot 
the intensity of emitted atoms from the lower band as a 
function of the decay rate of the optical resonator (kq) 
and the atom out-coupling rate (k^) for fixed atom-field 
coupling 77. The system is assumed to reach a stationary 
state due to atomic replenishing as described in more 
detail in the next section. We clearly see that a high- 
Q optical cavity enhances the atom emission almost as 
much as a low atomic loss rate. This is directly related to 
the dominance of stimulated pairwise atom-photon gen- 
eration. The corresponding intensity correlation function 
(pair correlation) 52(0) shown in Fig. [2Ib) clearly illus- 
trates the transition from thermal bunching 92 (0) = 2 to 
a Poissonian value 32(0) ~ 1 as for a laser crossing thresh- 
old. Atomic anti-bunching seems possible even without 
any direct interaction between the atoms. The stimu- 
lated enhancement of the twin gain via atoms and pho- 
tons is shown in Fig. where we depict the stationary 
upper-atomic-state population. Due to this enhancement 
the upper level gets depleted below one atom on average, 
and the atoms are preferentially deposited in the atom 
laser mode before they leak out. 

While the above numerical results confirm the basic 
idea, these calculations do not come near to experimen- 
tal realistic atom and photon numbers. Hence in the 
following we will develop further analytic approximative 
descriptions to reach more interesting and realistic pa- 
rameter regimes. 



III. INCOHERENT INJECTION OF ATOMS 

We first assume that the injection of atoms is realized 
by a pre-cooled atomic beam injected into the trap. This 




(c) 

FIG. 2: (a) Emission rate of condensed atoms from the lower 
vibrational band as a function of cavity linewidth Ka and atom 
out-coupling rate Hb for /iout ~ 0.9/iin ~ Sij/S and Kcb = rj/SO. 

(b) Intensity correlation function (72(0) of the condensed 
atoms from the lower vibrational band as a function of cavity 
linewidth Ka and atom out-coupling rate Kb for parameters as 
above. 

(c) Atomic occupation number of the upper (source) vibra- 
tional band as a function of cavity linewidth Ka and atom 
out-coupling rate Hb for parameters as above. 



will incoherently add and remove occupations from the 
atomic trap levels. Since the trap acts as a resonator cav- 
ity for matter waves, we model this by incoherent occupa- 
tion and depletion of the individual levels. In principle, 
we have four individual rates, altogether. However, as 
the ground mode is not intended to be directly pumped, 
we restrict the injection of atoms to the upper atomic 
(source) mode c. We further allow atoms to be lost from 
mode c due to heating or re-thermalization processes. 
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which might be included into the loss rate for the upper 
level. Identifying /Xi„ and /iout as the atomic incoming 
and outgoing rates, respectively, a corresponding model 
term in the master equation reads 



>CpumpP = Mill ( [c^ P c] + [cV, c] ) 
+ fJ-ont 



(6) 



Let us now have a look at the dynamics induced by 
this model. When there is at most one atom at a time 
present in the source mode, we can search for a solution 
in the reduced basis \n,m,k) representing n photons in 
the cavity mode, m atoms in the ground level, and fc = 
or 1 atoms in the source level. The joint probability of 
having n photons and m atoms is conveniently written 
down as 



with the diagonal density-matrix elements 
Pn.Tn,k = {n,m,k\p\n,m,k). 



(7) 



(8) 



The evolution of the joint atom-photon number statis- 
tics Pn,m{t) is hence determined by a set of coupled dif- 
ferential equations for the probabilities Pn,m,o and Pn,m,i 
which can be directly calculated from the master equa- 
tion. They read 



dt 



Pn,mfi = 'nV^Vn,ni 

+ 2Ka [{n + l)Pn+l,m,0 ~ npn^rnfl] 
+ 2Kb [{ni + l)p„,m+1.0 - mpn.nifl] 



(9a) 



■^Pn.mS = (n + + 1) 14+1, m+1 

+ 2Ka [{n + 1) 

7^71,771+1-1 ^ Pn,m,l\ 

+ 2n 

in Pn,m,0 2/iout Pri,m,li 



(9b) 



dt 



Mout )Vn .rn 



' Pn,'m,0^ 



+ 2Kay/n{n + l)K+i,m - (2n - 1)kqK 
-I- 2Kby^ m{m + l)K,m+i - {2m - l)KbVn 

where the off-diagonal elements are defined as 

Vn,rn ^ {u- 1,171- l,l\p\n,m,0) 

+ {n,m,0\p\n — 1, m — 1,1). 



(9c) 



(10) 



Note that in Eqs. ([9]) we have dropped terms involv- 
ing the dynamics of higher-order probabilities Pn,m,i>i 
to guarantee that the feed of atoms into the source mode 
has no direct influence on the atom-photon distribution 
Pn,m{i)- (In fact, this procedure is equivalent to replac- 
ing p„^m_2 by its steady-state value [lB|.) Mathematically, 



the evolution equation for Pn,m due to atomic injection 
must not depend on either of the parameters /ijn, /iout 
or, equivalently. 



I pump 



0. 



(11) 



In total, the evolution equation for the atom-photon 
statistics Pn.m obtained by adding Eqs. (Pa|) and (|9b[) 
reads 



V 



\/{n+ 1)(to + 1) T4+i,m+i 



2Ka [{n - 
2Kb [{m - 



l)-Pri+l,m - 
- l)-fri,m+l 



^ Pn,m\ 
7TL Pn.ra\ 



(12) 



We now want to derive a rate equation for Pn,m in 
the limit where the non-diagonal elements Vn^m can be 
eliminated. Whereas the relaxation time of Pn,m is on 
the order of k~^, k^^, Eqs. (O evolve on a timescale 
which is influenced by the atomic pumping dynamics, 
too. When /Lti„,/Xout ^ «^a,bi the source mode reaches its 
stationary state much faster than the laser modes after 
the creation of an atom-photon pair. The time scales 
are then significantly separated and one can adiabatically 
eliminate the fast evolving terms. 

The procedure is as follows. First, neglecting all terms 
proportional to Ka and Kb, the steady-state solution of 
Eqs. ^ is obtained by equating the time derivatives with 
zero. It is then with help of Eq. ^ expressed in terms 
of the slowly-varying probabilities Pn,m which yields 



(13) 



In this equation we have defined the emission and ab- 
sorption rates respectively as 



Br, 



nmR 
1 + nrriR ' 
nmR 



Pout 



1 + nmR 



(14a) 
(14b) 



and R = 2?]'^ / {pi^ + Pout)^- Finally, inserting Eqs. 113] 
and (|14p into Eq. we get for the rate equation 



dt 



Pn,m ^n-|-l,m-|-l-^7i-|-l,m-|-l -^^^-)-l,m+l^^^,^^ 

+ 2Ka [{n + l)P„+i,m ~nPn>m] 
+ 2Kb [{m + l)Pn,m+l ~ m 



(15) 



The interpretation of the individual terms is rather sim- 
ple. While the last two lines originate respectively from 
photon loss and atomic out-coupling, the first two lines 
describe absorption (terms Sn,m) and creation (terms 
~ An^m) of laser atom-photon pairs. This equation is 
in fact very analogous to the rate equation found for a 
common single-mode laser showing gain and saturation. 
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(n'=m') =Yjnni n'^m'Pn 



FIG. 3: (a) Photon and atom mean occupation numbers 
(n) = (m) in the steady state as functions of the atomic pump- 
ing rate /iin- They coincide since the decay rates equal and 
Ka = Kb = r?/20. The solid line and the (red) dotted line rep- 
resent, respectively, the numerical solution of the rate equa- 
tion (|15p and the approximate analytical solution from Eqs. 
(|17[) showing good agreement sufficiently far above threshold. 

(b) Single-mode number distribution P„ (crosses) for the same 
parameters as in (a) for fiin = 50k and, for comparison, a Pois- 
son distribution with equal mean (dotted line). 

(c) Single-mode correlation function g^^\0) (solid line) and 
inter-mode correlation function g^'"^ (dashed line) as func- 
tions of the pumping rate /iin for the same parameters as in 
(a). The dotted lines represent the corresponding analytical 
results. 

(d) Joint excitation probability P„,m for the same parameters 
as in (b). The quenching along the broken line is significant 
for enlarged inter-mode correlations. 



As a central result, however, we find here joint atom- 
photon gain and saturation in the emission and absorp- 
tion rates [Eqs. (HH)]. The origin for this twin Bosc- 
enhancement which has aheady been anticipated |23l . l3l| 
is cavity-assisted cooling of atoms as reahzed by Eq. 
The simultaneous creation of ground-level-atom-photon 
pairs is further responsible for the symmetry between n 
and m in Eq. (jlSp . For equal decay rates Ka — we 
have perfect symmetry Pn,m — Pm,n and both lasers fea- 
ture identical properties. 

In the following we are interested in the stationary so- 
lution of the system which can be obtained by taking the 
time derivative in Eq. (|15p to zero. Before solving the re- 
sulting set of equations, let us first try to obtain analyti- 
cally a solution sufficiently far above threshold where we 
consider the probabilities Pn,m different from zero only 
for large n and m. We are then allowed to replace An^m 
and Bn,m in Eq- ([T5)) respectively by /Zin and /Zout and 
get the following dynamical equations for the moments 



[m) 



2Ka{n) -I 
2Kb(TO) - 



Mill - /^out, (16a) 
-Miii-/^out, (16b) 

f 2 (^in - A*out + Ka) (n) 



+ Miii+Mout, (16c) 
(to2^ ^ _ 4^^(771^) -I- 2 (^in - Hout + Kb) (m) 

+ Miii + Mout, (16d) 
(nm) — — 2 (na + Kb) {nm) 

+ i^J■m ~ Mout) i{n) + (m)) + + ^J,out■ (16e) 

The threshold condition is here //in > /iout- Note that 
Eqs. (|16p are only valid when — Mout ^ 2 min (kq, Kb) 
and should then reproduce the semiclassical results. In 
the stationary state, the mean photon and atom numbers 
are 



/^in Mout 

2Ka 
/^iii /^out 
2Kb 



(17a) 
(17b) 



As expected they coincide when Ka = Kb- In the remain- 
ing of this section we will discuss this situation. We will 
therefore write k instead of Ka and Kb and refer to both 
modes simply by n, too. 

Figure[3l^a) shows the mean occupation number (n) for 
K — 0.072i? = 0.1/iout (dotted line) as a function of the 
pumping strength /iin together with the numerical result 
obtained directly from Eq. ([15]) (solid line). Sufhciently 
far above threshold we get perfect agreement. More no- 
tably, however, both modes are in a state very close to 
a coherent state. This is illustrated in Fig. [3Jb) where 
the single-mode statistics P„ = J2m Pn,m is plotted (+) 
compared to a corresponding Poisson distribution with 
equal mean value (dotted line) for fi-m/ k = 50. 

It is interesting to ask whether or not our system shows 
nonclassical features such as photon anti-bunching or 
noise reduction. As a measure for the number statis- 
tics we compute the second-order coherence function at 
r = 0, 



(18) 



where (An^) denotes the photon-number dispersion. The 
steady-state solution of Eqs. (fTB|) yields 



9'^'\0)^1 + 



Z^out 1 

~2i^'M' 



(19) 



which is represented by the (blue) dotted line in Fig. 
Wis). Again, we have plotted the numerical result from 
the rate equation for comparison (solid line) and find 
agreement for large /iin. There the correlation function 
approaches unity, the value for a coherent state. How- 
ever, (/(^^ (0) > 1 even far above threshold, which excludes 
photon anti-bunching being reached. 
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Nevertheless, there are nonclassical effects apparent. 
Let us first have a look at the inter-mode correlation 
function 



gn,ra 



(n) (m) 



(20) 



which we can derive from the steady-state solution of 
Eq. (fT5| . We have depicted the result as the dashed line 
in Fig.jSjc). It is well known [39|], that the inter-mode cor- 
relations exceed the correlations of photons of the same 
beam in two-photon lasers. Mathematically, this implies 
violation of a corresponding Cauchy-Schwartz inequality 
which reads in our situation (due to the symmetry) 



(21) 



lies 



Note, that here we recover the same feature: g' 
always above 5^^^(0) in Fig.jSlJc). 

By use of the steady-state solution of Eqs. we can 
estimate that 



g".™ ^ g(2)(0) 



Mout 



2k 



(22) 



which indeed proves this violation, at least sufRciently 
far above threshold. The (red) dotted line in Fig. [3^c) 
corresponds to the solution of g"'™ as derived in Eq. ([22|) . 

Another interesting feature in twin-photon lasers is 
noise reduction of one quadrature of the field below the 
classical level [1^. Perfect squeezing in the stationary 
output-beam intensity difference was found [4l| for which 
simultaneous creation of the twin photons was identified 
to be the origin. 

Let us see whether we can find here any similarities. 
Figure [3Kd) gives a contour plot of the joint probabil- 
ity distribution Pn,rn in the steady state indicating the 
afore mentioned strong inter-mode correlations. Due to 
the elliptic shapes of the curves we anticipate a reduced 
width of the intensity difference, i.e., along the broken 
line. From the rate equation (fTS]) we deduce that 



dt 



(n — m) = — 2K(n — m) 



(23a) 



dt 



{{n - mf) = -2k [{{n - mf) - {n + m)] (23b) 



such that the width in the intensity difference is given by 

((n-m)2)_(n-m)2 = ili±^, (24) 

both above and below the oscillation threshold. On the 
other hand, classical theory predicts the inequality 

((n - m)2) - (rt - m)2 > (n + m) (25) 

and we thus have noise reduction of 50% below the clas- 
sical level. 



IV. ATOMIC PUMPING BY TUNNELING 

Based on atomic tunneling between lattice sites, we 
will present in what follows a model for atomic injec- 
tion which provides a proper and consistent description 
of pumping of the system. We allow here only atoms to 
move within the excited band since, in general, the rates 
for atomic tunneling involving the lower band are usu- 
ally much smaller. Let denote the creation operator 
for upper-band atoms within the cavity and the operator 
r represent the available lattice sites outside the cavity 
which are coupled to the inside. A corresponding tunnel- 
ing Hamiltonian is 



Ht = ic'T 



(26) 



where F, besides reducing the number of atoms outside 
by one, also includes the interaction strength. Let us 
emphasize here that Eq. ((26)) can in principle also be 
viewed as an interaction changing the atomic state from 
an untrapped state to a trapped state. Atoms in the un- 
trapped state would in this situation populate the mode 
represented by F and be transferred to the source mode c 
by a proper interaction. This interaction could, e.g., be a 
Raman transition applied on untrapped atoms which are 
already present in the cavity. However, to get a suitable 
description of pumping, we will include Ht into the mas- 
ter equation [Eq. ([3])] and then define a system density 
operator by taking the trace over the levels represented 
by F. 

First, we have to extend our three- mode system Hilbert 
space, which we may denote as Hsys — span{|n)}„=i...oo, 
by the space Hont spanned by the basis {\J)}j=i...oo, 
where |J) denotes a superposition state containing J 
atoms. A general state can be written as \n, J) = 
\n) (g) I J), where |n) G Hsys and |J) G Hout- For the 
total density operator we use the notation 

Ptot = ^ XnrnJK\n,J){ni,K\ (27) 
nrnJK 

and define our system density operator as the trace 

p = ^(/|ptot|/) =^A„„|n)(m| (28) 

/ nm 

with coefficients given by the sum over the diagonal ele- 
ments, A„„i = J^j'^nmjj- After incorporating the tun- 
neling Hamiltonian into the master equation and tak- 
ing the trace over Hout, the pumping term reads 



(29) 



where a = J2 XnmjK {K\r\J)\n) {m\. 

In the remainder of the paper we will consider the 
scenario where the number of outside modes is rather 
large and, further, these modes are in a state which is 
marginally affected by atom transfer to or from the in- 
tracavity system. (For example, this holds true for a co- 
herent state or for a state with many sufficiently occupied 
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levels.) It is then justified to approximate a ~ /zp, where 
/i indicates the effective pumping strength, so that from 
Eq. ((29|) we can read off an effective pump Hamiltonian 
of the form 



Hp = ifj, (c^ — c) 



(30) 



Note that without loss of generality we have assumed here 
fi to be real thus having zero phase. It is of course an 
artifact of the model that the atomic field exhibits here 
a phase. However, Eq. ([26]) conserves the atom number 
such that the phase has no observable consequences. Hp 
thus provides a self-consistent description of continuous 
atomic occupation of the source mode originating from 
tunneling. 

In addition to pair generation [Eq. ([2])] and atom injec- 
tion [Eq. ((30|) ] let us also include atomic heating into the 
dynamics. Heating should regard both the atomic laser 
mode — atoms are heated into the source mode at rate Kbc 
there — and the source mode, which we model simply by 
a depletion channel (rate Kc). The full master equation 
reads then 

P = - i[Hint + Hp,p] 

+ Kt{[b,pb^ + [bp,b^]) 

+ Kbc{[c^b,pb''c\ + [ct6p,6tc]) 

+ «:e([c,pct] + [cp,ct]). (31) 



We proceed by converting the master equation ((3T|) 
into a Fokker-Planck equation (FPE) for the generalized 
P-representation of the density operator. The general- 
ized P-representation is defined as 



p^ [ dQdatp(Q,Qt) l°^H°^^ 
J (ot^ \a) 



(32) 



where a = (a,/3, 7) and — (a^,/3^,7^) are inde- 
pendent c omp lex variables. Applying standard operator 
identities [42|, the FPE for the P-representation reads 

d d 



d_ 

"dp 



[-Ktf3 - Kbcf3 (7^7 + 1) + '7 a^l] Pia, a^) 



+ „,-, ■>] 7 -P("- "^) + aoQ i^l^bcP 7) P{(^, Oi^) 



dadp 



92 



Kb, 



8/38^ 
;/3'f/?P(a,Q:t) -l-conj. 



(33) 



where the conjugate terms are given by the formal re- 
placement a <-> a^. The use of the P-representation en- 
ables us to derive normally ordered expectation values of 
quantum operators simply as the corresponding moments 
of P(q;, Q;t). However, we will not present here a solution 



to the FPE ((33)) but use an equivalent set of stochastic 
differential equations (SDE) for the variables a and . 
Therefore, we can identify quantum expectation values of 
any polynomial p composed of the operators a = (a, 6, c) 
with the mean of our stochastic variables as 



(:p(a,at):) = (p {a,a^))^ 



(34) 



where the double dots indicate normal ordering. As an 
indication for normal order we will later express the quan- 
tities of interest in terms of the intensities Jq. = a^o;. 
(Note that the bold variables are used here as a formal 
placeholder for one of the variables a, /?, and 7 at a 
time.) This causes normally ordered expectation values 
to be calculated simply as 



(35) 



where fia denotes the corresponding quantum number 
operator. 

The SDE have the same drift and diffusion terms as 
the FPE §^ and read 

a = —KaU + r](3'^^ + S,a, (36a) 

13 = -Kbp - Kbcl3 (7^7 + 1) + + ip, (36b) 

7 = - (kc - K6c/5^/3) 7 - '7a/3 + M + (36c) 

where the ^q, are (S-correlated noise terms with zero mean 
and the following non-zero diffusion coefficients. 



'Dad ^vil), 

^'at/Jt = 77(7^), 
15/37 = -KbciPl), 
T^I3ij^ = -KbciP^l^) 
V^,^ = KbciP^lS). 



(37a) 
(37b) 
(37c) 
(37d) 
(37e) 



Before we proceed, let us first investigate the semiclas- 
sical behavior of the system. It is gained by ignoring the 
stochastic nature of the system variables and equations, 
i.e., by taking = a* and = in Eqs. 

The steady-state solution (a) = ag can be obtained by 
equating the drift terms in Eqs. p6p with zero. As men- 
tioned, we chose for convenience the pumping strength 
p to be real so that the semiclassical solution olq has no 
phase. The result shows a threshold appearing at 



Ka{Kb + Kbc) 

Mth = Kc\ 5 

\ V KaKbc 



(38) 



for > KaKbc, below which the mean field amplitudes 

are aQ — Pq — and 70 = p/kc, while above we have 



ao 



Tj KaKbc 



V\/KciKb + Kbc) , . 

7, ve-1, (39a) 

(39b) 

70 = — , (39c) 



„ / ^'CL I 

Po = ve - 1 

r\ KaKbc 
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where e = fi/fith is the pumping parameter. The mean 
occupation numbers are then = ckq, respectively. 

The semiclassical solution (a) may provide an approx- 
imation to the solution of the FPE ([55]) for small fluctua- 
tions. As an example, for lasers this is fulfillcd far above 
threshold [i^l. Nevertheless, it does not account for cor- 
relations of any kind and, in general, high-order expecta- 
tion values of a disagree with the semiclassical counter- 
part. One is therefore obliged to include the fluctuations 
into the dynamics in order to calculate for instance the 
number statistics or output spectra of the laser. 

In what follows, we will again take the noise terms 
^Q- into account and derive equations for the dynamics 
of the fluctuations around the steady-state solutions chq. 
These quantities turn out to be directly related to the 
noise spectra from which we can calculate higher-order 
moments in an appropriate way. 



V. DYNAMICS OF FLUCTUATIONS 

To account properly for the fluctuations, we now ex- 
press the dynamic variables as a sum of the steady-state 
values plus small fluctuations. 



a = cxq + Sa, 



(40) 



(and analogously for a^) and linearize the equations of 
motion around their steady-state values. Note that it is 
important to be in a parameter regime where the fluctu- 
ations of the dynamical variables are much smaller than 
their steady-state values, in order to justify this assump- 
tion. We will later see that this criterion can be fulfilled 
even not very far above the threshold. 

After substitution of Eq. ^ into the SDK the 
linearized equations for the fluctuations Sa and Scy.^ are 



d 



6a 



= M 



Sex 



where the matrix M is composed of four blocks. 



M 



" A 


B ' 


B 


A_ 



(41) 



(42) 



with 



A = 



and 







B 



77/3o 



(43) 



7770 
Kbc/3o7o 








(44) 



The solution of Eq. (|^T|) is best found by a Fourier 
transform method. Let us define the variables v = 



{Sot, 5a^ , S^a, Cat) in frequency space as 

v{lu)^^J dre— t;(T) 



(45) 



and note that we will for convenience use the same sym- 
bol both in time space and in frequency space while, 
though, the argument should make the distinction ob- 
vious. It follows then from this definition that the noise 
terms ■Cct(w) obey the correlations 



(46) 



with the same diffusion coefficients as the untransformed 
variables given in Eqs. (j37p . Applying the Fourier trans- 
form on Eq. (|¥T1) yields in the steady state 



6cy.{uj) 
6a'' {lu) 



(47) 



which is an algebraic problem that can be solved exactly. 
The solution in time space is then — of course — governed 
by the inverse Fourier transform. 

Using Eqs. and (|T7)) . we define the spectral func- 
tions AafB^Lu) and Qap{^) written in terms of the quadra- 
ture component 5xa = 5a + 6a^ as 



{6a\u)5p{Lo')) ^ Ao,p[uj)5{Lo + u'), 



{5xo,{uj)6xi3{uj')) 



iaf3 



{lu)6{uj + uj'). 



(48a) 
(48b) 



They are related to the spectra of, respectively, the am- 
plitude and quadrature fluctuations by just adding the 
corresponding shot-noise contributions. We will come 
back to this point later when we discuss number statistics 
and squeezing. 

According to Eq. ([^0]) the mean populations ria = (la) 
are given by 



a'^a) = an 



(Sa^Sa), 



(49) 



where we have used the fact that {6a) ~ {Sa'') = 0. The 
correction to the semiclassical solution is obtained by the 
integral over the amplitude fluctuation spectrum. 



{Sa^'Sa) — — 
^ ' 27r 



dw Aq,oi(w). 



(50) 



The result for the photon laser intensity (/„) is shown as 
the dash-dotted (blue) fine in Fig. [4{a) as a function of 
the pumping parameter e. The remaining parameters are 
(r?, Kbc, Kc) = (4, 0, 2)kq. Therefore the solutions for both 
modes a and h coincide for equal decay rates = k^. 
The mean number follows closely the linear semiclassical 
result which we have indicated by the crosses. 

However, care must be taken for [i k, /ij^ since the fluc- 
tuations around the steady-state may be there too large 
to justify the linearization assumption. For comparison 
we have also plotted the quantity {Sa'Soi) as the dotted 
(red) line which turns out to be much less than one even 
not so far above the threshold. 
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A. Number statistics for equal decay rates 

Laser light is distinguished by having smaU intensity 
fluctuations in the sense that the standard deviation of 
the intensity is small compared to its mean. The fluctu- 
ations are usually investigated by means of the variance 
of the photon number. It is well known that the light 
emitted from a laser operating far above the oscillation 
threshold exhibits the statistics of independent particles, 
which corresponds to a Poisson distribution for the num- 
ber of photons in the laser mode. 

Let us write down the photon-number variance in nor- 
mal order which reads 



{Anl 



(:Ahi 



(51) 



The first term on the right-hand side of this equation 
originates from the commutator of the bosonic opera- 
tors a and and reflects the quantum character of the 
fleld. In fact, it corresponds to the contribution to the 
number variance of an ensemble of independent parti- 
cles. Commonly called the shot noise, it is in principle 
inherent to any state of the electromagnetic field due to 
the discreteness of the photon-number distribution. For 
Poissonian number statistics it constitutes the sole contri- 
bution to the number variance, and field states for which 
the photon-number variance is less than the shot-noise 
term are necessarily nonclassical '43] . 

As a consequence of using the P-function represen- 
tation, which is mathematically manifested in equality 
([55)1 , the second term on the right-hand side of Eq. ([5T|) 
is equal to the c-number intensity dispersion. 



■■Ahl 



where 



Hi 



(52) 



(53) 



This term exhibits the characteristic of the P-function 
playing the role of a quasi-classical distribution function. 
We can hence associate the deviation from the Poisson 
distribution with the intensity dispersion so that negative 
and positive (AJ^) indicate sub- and super-Poissonian 
statistics, respectively. 

Intensity-correlation experiments reveal the cumula- 
tion tendency of photons of a light beam. Usually, one 
is interested in the joint probability of detecting a pho- 
ton respectively at times t and t + t. In the stationary 
state, the probability is proportional to the second-order 
coherence function 



i{t)n{t - 



(54) 



Note, that one would have g^^^r) = 1 for independent 
particles while g'^-'(O) < g'-^-'ir) is significant for photon 
anti-bunching. 

On the other hand, there exists a relation between the 
degree of second-order coherence and the number statis- 
tics. Obviously, the numerator in Eq. (j54p is equal to 




(fot<5a> 



/ ■ 



.9f'{0) 



(a) 





11 16 0.1 1 10 100 



(b) 



FIG. 4: (a) Single-mode correlation function (0) (solid 
line) and inter-mode correlation function g^'^ (dashed line) 
where the variable is the pumping strength \i. In addition we 
have plotted the mean intensity (Jo) (blue dash-dotted line) 
together with the corresponding semiclassical result (-I-), and 
the quantity (5q:^(5q) (red dotted line). 

(b) Dependence of the correlation functions (j^-* (0) (solid line) 
and g°''^ (dashed line) on the out-coupling rates Ka = Kb- The 
pumping strength /x is here adjusted to yield mean intensities 
of around one as indicated by the (red) dotted line. 



the normally-ordered expectation value of if r = 0. 
It is thus possible to re-express the coherence function 
for r = in terms of the number variance. As before we 
use here the intensity variance (to account for the normal 
order) and get 



5i^H0)-l + 



(Ai| 



(55) 



So far no approximation has been made. The intensity 
dispersion in Eq. (j55p involves fourth-order correlations 
of the field variables a. Nevertheless, in a regime of small 
fluctuations, one can compute this quantity in terms of 
the amplitude quadrature variance which requires only 
second-order correlations. Analogous to Eq. (|40)) we 
therefore write for the intensity = (/„) + SI^ so that 
the variance in the intensity is equal to (((5/q,)^). We fur- 
ther assume that the classical relation between intensity 
and amplitude fluctuations 



Sir. 



cxq (^a^ + (5q;) = cxoSxc 



(56) 



holds true. Note that this approximation is justified as 
long as (Sot^Scy.)'^ <C ot^{Sx'^). We get in this case for the 



intensity dispersion 

where, written in terms of the spectral functions 



(57) 



(58) 



When |(^a;^)| ^ 1, the intensity dispersion can be ne- 
glected compared to its mean. According to Eqs. ([?T|) 
and (|52p . the photon number variance is in this case 
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(An^) = (no,) which is typical of Poissonian number 
statistics. Equivalently, the second-order coherence func- 
tion at T = equals one, which follows from Eq. ([55|) . 
The P-function is then very close (or equal) to a S- 
distribution representing the density matrix of a coher- 
ent state. This behavior is demonstrated in Fig. \M^a,) 
sufficiently far above threshold. The solid line shows the 

coherence function ga (0). We have chosen here a set 
of parameters such that the solutions for modes a and 
b coincide, i.e., Ha — Hb — k.cI'^ = '7/4. At and slightly 
above threshold, the system is in a state close to a ther- 
mal state and (0) < 2. With increasing pumping 
strength the laser modes start to become populated. It 
is significant for the stimulated enhancement of the de- 
sired processes that the coherence function drops at the 
same time — quite sharply — to its asymptotic value. For 
the parameters chosen here it is unity which corresponds 
to a coherent state. 

In addition we have also plotted the inter-mode corre- 
lation function, which is defined as 
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(ata)(/?t/3)^ 



(59) 



as the broken line. As a manifestation of the pair creation 
we find large correlations between atoms and photons. 
Actually, the curve for g°-^ lies always above the curve 
that belongs to the coherence functions, which implies 
that the mutual correlations between atoms and photons 
are much larger than the correlations between the parti- 
cles of one and the same beam are. 

Let us mention at this point that Eq. (|56p relates num- 
ber squeezing to amplitude quadrature squeezing. Am- 
plitude and number squeezing are not equivalent. As an 
example, amplitude squeezed states can have both super- 
and sub-Poissonian number statistics for fixed photon 
number, depending solely on the squeezing parameter 

m. 

Equation (|58p is proportional to the normally ordered 
quantum expectation value of the variance of the field 
amplitude quadrature X = a -I- at , namely 



(5a;i) = (:AX^ 



(60) 



Amplitude quadrature squeezing, where (: AX^ :) < 0, 
is therefore only possible when the corresponding P- 
function is not positive definite. As a consequence, the 
intensity dispersion (|57[) is negative and the modes ex- 
hibit sub-Poissonian statistics. Mathematically, the co- 
herence function (0) < 1. 

This situation can be achieved by reducing the decay 
rates as is demonstrated in Fig. I^b). It shows the coher- 

ence function g)i (0) (solid line) and the correlation func- 
tion g""'^ (broken line) for ki, — Ka- As already mentioned 
in the discussion of Fig. O reduction of the out-coupling 
rates induces stimulated enhancement of the twin gain 
with large atom-laser emission rate for a cavity with high 
quality. In principle, we find here the same results. How- 
ever, here we would like to keep the mode occupations 



constant for a more refined comparison. Therefore the 
pumping strength fi is adjusted here to achieve mean 
intensities (la) of the order of one by use of the semi- 
classical result. In the figure, (/„) is indicated by the 
dotted line. The other parameters are rj — IOkc and 
Khc = so that the results for modes a and b coin- 
cide. For large Ka, the second-order correlation function 
starts from two where the mode is in a thermal state. It 
decreases when the linewidth is reduced and eventually 
drops below unity. We have checked that at the same 
time the threshold pumping strength //th decreases lin- 
early with the cavity line width. A good cavity thus not 
only enhances the emission of ground-level atoms but also 
eases the requirements on the pumping strength (i.e., the 
atomic density) to reach degeneracy. 

Note finally that the inter-mode correlation g""'' is 
again always larger than g^^^ (0) and reaches an asymp- 
totic value of 1. 



B. Output-field spectra 

So far we have focused on the internal properties of our 
system and were able, for instance, to derive the mean 
number and the statistics of the laser modes. In the scope 
of any experiment, however, one is usually interested in 
accessible quantities such as the output fields. We will 
therefore discuss now the properties of the output beams, 
particularly the fluctuation spectra of the output fields 
and of the intensity difference. 

The result (j56|) which applies for instance for an laser 
operating far above the threshold not only connects 
quadrature and number squeezing. Moreover, it gives 
a relation between the number statistics of the inter- 
nal fields and the spectra of the output fields which are 
accessible in experiments. The fluctuation spectrum in 
the field is given by the Wiener-Khinchine theorem [i^ j. 
We will use here the input-output formalism developed 
by Collet and Gardiner [i^ to express the output-field 
operators 6out(i) in terms of the system variables. The 
spectrum of fiuctuations corresponding to the quadrature 



b'^^^ can then be written as 



Vaiuj) - 1 + 2k„ / dTe'^^6x^{T)SXo.{0)). 



(61) 



The first term on the right-hand side of this expression 
originates from the commutator of the operators bout and 
6^^^ when putting them into normal order and corre- 
sponds to the shot noise. Equation (pT|) may also be 
expressed in terms of the spectral functions Eqs. 
which yields 



Vai^j) = 1 + 2Ka(3aa(^)- 



(62) 



As we have seen previously, sub-Poissonian photon statis- 
tics is a consequence of amplitude-quadrature squeezing. 
According to Eqs. and (p^ . one must have in this 
case Va{i-^) < 1, at least in some frequency range. This 
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(a) 



10 20 



(b) 



FIG. 5: (a) Fluctuation spectrum Va(cj) for different values 
of the pumping parameter e. The areas below the curves 
are directly correlated with the number statistics (see text). 
Sufficiently far above the threshold (solid line) the system 
exhibits amplitude noise reduction and sub-Poisson statistics, 
(b) Spectrum of the intensity-difference fluctuations in the 
output beams, ^(aj), for = (dashed line) and = 
(sold line). The inset shows the contribution at zero frequency 
S'(O) as a function of the heating rate Khc- 



type of noise compression is illustrated in Fig.[5l^a) show- 
ing the fluctuation spectrum T4(ijj) for different values of 
the pumping parameter e. The other parameters are the 
same as in Fig.UJJb) with Ka = Kc/5. Close to the thresh- 
old, the spectrum 14 (w) lies above one as is indicated by 
the dotted line for which e = 1.5. With increasing pump, 
noise reduction in the low-frequency range occurs until 
eventually the contributions of noise attenuation and ac- 
centuation cancel. This is illustrated by the gray areas in 
the plot for which e = 2.7 and (/'^•'(O) = 1 (dashed line). 
Far above the threshold the statistics is sub-Poissonian, 
the solid line in the figure corresponds to e = 20 with 
5(2) (0) = 0.85. 

As a result of simultaneous atom-photon pair genera- 
tion, the atom and photon numbers in the lasing modes 
are strongly correlated. They are, actually, exactly equal 
so long as either an atom or a photon is removed from 
the system. The variance of the occupation-number dif- 
ference vanishes therefore in a system without losses. 

For a steady operation including reloading and deple- 
tion of atoms we already found in Sec. IIIII a reduction of 
the variance in the number difference. Provided that all 
atoms and photons are detected, this feature should be 
recovered for the output intensities, as is the case for a 
parametric amplifier. In parametric amplification which 
our system is very closely related to, perfect noise reduc- 
tion at zero frequency was found in the output spectrum 
of intensity-difference fluctuations. 

Let us see what this means here. The fluctuation spec- 
trum in the output intensity-difference is given by 



Auj e" 



(A(r)A(O)) - (A)^ 



(63) 



where A 



-''out "out 



&Q^j6out- As before we can express 
the spectrum in terms of the internal variables and get, 
after putting the right-hand side of Eq. (j63p into normal 



order, 

S(ijj) = 2Ka{Ia) + 2Kb{Ib) + (2Kaao )^(3aa (w) 

+ {2KbPofQ0[3{uj) - 8KaKbaof3oQa0iuj). (64) 

Here, the first two terms are the shot noise. Figure 
[5jb) displays the normalized fluctuation spectrum (j64p 
divided by the shot-noise contribution. When the heat- 
ing rate Hbc = vanishes (dashed line), we find perfect 
noise reduction at resonance. Even for unequal output- 
coupling rates Kq 7^ Kf,, the zero- frequency contribution 
5(0) vanishes and the system still realizes a perfectly 
quantum correlated pair source. 

Nevertheless, atomic scattering due to coUisional heat- 
ing or photon emission into modes other than the cavity 
mode can cause deviations in the properties of the out- 
put fluxes. For an illustration we have chosen the atomic 
heating rate Kbc as high as the out-coupling rates Kq = Kb 
for the solid line in the plot. While there is still noise 
reduction in the low-frequency range, the noise may be 
slightly increased for higher frequencies (not shown in 
the plot). The contribution at resonance is shown in the 
inset as a function of the heating rate Kbc- 



VI. CONCLUSIONS 

We have shown that by help of an optical cavity one 
can implement the final cooling stage of an atom laser via 
Raman sideband gain. As the corresponding anti-Stokes 
transitions are enhanced by stimulated emission through 
the photons and atoms present, the setup should allow 
for the realization of an atom-photon twin laser. Gener- 
ation of cavity photons is always accompanied by atomic 
transitions to the atom laser mode, such that photon 
emission through the mirror acts as a diode for atom ac- 
cumulation in a single state. The more photons present 
the more atoms are transferred to the lower state on av- 
erage. Hence this device produces two output beams con- 
sisting of photons and atoms, respectively, with genuine 
properties. 

Compared to conventional sideband cooling, the use 
of a cavity reduces the requirements on the atomic den- 
sity (threshold for degeneracy) and avoids the problems 
of photon spontaneous emission and reabsorption as well 
as atom-atom interactions which are detrimental to the 
atom-laser phase [l8| . When quantum degeneracy is 
reached, an atomic beam is provided by a proper out- 
coupling mechanism. 

Using a simple model for incoherent atomic pumping 
we were able to derive from the master equation a rate 
equation for the joint atom-photon number distribution 
very analogous to the usual laser equations, but even 
more reminiscent of an optical parametric oscillator op- 
erated in the strong coupling limit [46|. Our model al- 
lowed to explicitly calculate absorption and creation rates 
of atom-photon pairs. Analogous to a laser the above- 
threshold population of both laser modes yields a close 
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Poissonian number statistics in the steady state with even 
some anti-bunching possible for the atoms. The corre- 
lations between particles and photons are significantly 
larger and even show squeezing in the number difference. 
In analogy to a twin photon source we thus could get a 
heralded single atom beam even if threshold cannot be 
fully reached. 

A central technical problem to solve is to achieve a suf- 
ficiently fast supply of new atoms. Any realistic descrip- 
tion of this of course would require at least 2D modeling 
of the atomic motion. While this seems straightforward 
in principle, such models would be tedious to solve in 
practice and require large simulation efforts. However, 
we do not expect qualitative changes in the results over- 
all. 

Let us finally mention that intracavity Raman gain 



has been already observed and used for cooling in an 
experimental setup by Vuletic and co-workers in Cesium 
[43] • Here the initial atomic temperature was quite high 
and there was no external lattice or trap present, but the 
results should give hope for current improved efforts 
starting from much colder sources. 
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